Cosmological models described by a mixture of van der Waals fluid and dark energy 
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The Universe is modeled as a binary mixture whose constituents are described by a van der Waals 
fluid and by a dark energy density. The dark energy density is considered either as the quintessence 
or as the Chaplygin gas. The irreversible processes concerning the energy transfer between the van 
der Waals fluid and the gravitational field are taken into account. This model can simulate: (a) an 
inflationary period where the acceleration grows exponentially and the van der Waals fluid behaves 
like an inflaton; (b) an inflationary period where the acceleration is positive but it decreases and 
tends to zero whereas the energy density of the van der Waals fluid decays; (c) a decelerated period 
which corresponds to a matter dominated period with a non-negative pressure; and (d) a present 
accelerated period where the dark energy density outweighs the energy density of the van der Waals 
fluid. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

One of the main objectives of the cosmological models 
is the description of the different phases of the Universe 
concerning the time evolution of its acceleration field. 
The first epoch represents a rapid expansion of the Uni- 
verse which is known as the inflationary period. Most 
of the theories that describe this period makes use of 
a scalar field which is related to a hypothetical parti- 
cle, the so-called inflaton (see e.g. the works pj and the 
references therein). The second period refers to a past 
decelerated epoch where the energy densities of the ra- 
diation and matter fields outweigh the scalar field. The 
present epoch is characterized by an accelerated Universe 
dominated by a dark non-baryonic matter and a dark 
energy density. The most common models for the dark 
energy density make use of the equations of state of the 
quintessence Q and of the Chaplygin gas Q • 

Another model for the Universe was proposed recently 
by Capozziello and co-workers Q who used the equation 
of state of a van der Waals fluid in order to analyze the 
accelerated behavior of the Universe. The advantage of 
this model is that it can describe the transition from a 
scalar field dominated epoch to a matter field dominated 
period without introducing scalar fields. 

The objective of the present work is to investigate a 
Universe described by the van der Waals equation of 
state. In section 2 we follow the works [4( and model 
the Universe as a non-dissipative one-component van der 
Waals fluid. We determine the time evolution of the ac- 
celeration, energy density and pressure fields and show 
that this model can simulate the two phases of the Uni- 
verse beginning with an accelerated period and going into 
a decelerated epoch. Although this model could describe 
the transition from an inflationary period to a matter 
dominated epoch it could not simulate the present ac- 
celerated period of the Universe. For this end we have 



modeled in section 3 the Universe as a binary mixture of 
a van der Waals fluid and dark energy. The dark energy 
is also regarded either as the quintessence or as the Chap- 
lygin gas. Furthermore, we have consider the irreversible 
processes related to the energy transfer between the van 
der Waals fluid and the gravitational field, since it is very 
questionable to get rid of the dissipative effects during the 
evolution of the Universe. Among other results, it was 
shown that this model can simulate: (a) an inflationary 
period where the acceleration grows exponentially and 
the van der Waals fluid behaves like an inflaton; (b) an 
inflationary period where the acceleration is positive but 
it decreases and tends to zero whereas the energy den- 
sity of the van der Waals fluid decays; (c) a decelerated 
period which corresponds to a matter dominated period 
with a non-negative pressure; and (d) a present acceler- 
ated period where the dark energy density outweighs the 
energy density of the van der Waals fluid. Units have 
been chosen so that c = 1 and 8nG/3 = 1. 



II. UNIVERSE AS A VAN DER WAALS FLUID 

In this section we consider a model for a spatially flat, 
homogeneous and isotropic Universe where the energy- 
momentum tensor of the sources of the gravitational field 
is described by a perfect fluid with a van der Waals equa- 
tion of state in the absence of dissipative processes. The 
energy-momentum tensor T^ v of the sources is given by 



(1) 



where (such that U^U^ = 1) is the four- velocity. The 



pressure of the van der Waals fluid p u 
energy density p w by (see, e.g. @) 



Pu 



&w w p u 



3p u 



is related to its 



(2) 
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In the above equation the pressure p w and the energy 
density p w are written in terms of dimcnsionlcss reduced 
variables and w w is a parameter connected with a reduced 



2 



200 - 



a 



-2110 




0.5 




FIG. 1: Acceleration vs time for Universe as a van der Waals 
fluid. 



FIG. 2: Energy density (straight line) and pressure (dashed 
line) vs time for Universe as a van der Waals fluid. 



temperature. Here w w will be identified with the coeffi- 
cient of the barotropic formula, since for small values of 
p w we have p w oc w w p w . 

By considering the Robertson- Walker metric and a co- 
moving frame one can get from the conservation law 
T^ v = for the energy-momentum tensor the following 
balance equation for the energy density 



3H(p w +p w ) = 0. 



(3) 



In the above equation H = a(t)/a(t) denotes the Hubble 
parameter, a(t) is the cosmic scale factor and the over-dot 
refers to a differentiation with respect to a dimensionless 
time t. 

The connection between the cosmic scale factor and 
the energy density is given by the Friedmann equation 



H 2 = p. 



(4) 



whose differentiation with respect to time leads to the 
time evolution of the cosmic scale factor 



H 2 
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3-H 2 
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(5) 



The solution of the second-order differential equation 
(J5J for a(t) can be found by specifying initial values for 
a(i) and H (t) at t = 0, for a given value of the parame- 
ter w w . Here we have chosen the values a(0) = 1 for the 
cosmic scale factor and H(Q) = 1 for the Hubble param- 
eter (by adjusting clocks) and in the following we shall 
comment on the possible values for the parameter w w . 

In the interval < w w < 1/2 the energy density grows 
with time, the pressure is always negative whereas the 
acceleration is always positive. One can conclude that 
this solution does not model properly the evolution of 
the Universe. 

The case where w w = 1/2 is the most interesting and 
in figures 1 and 2 we have plotted the acceleration, the 
energy density and the pressure as functions of the time. 
One can infer from figures 1 and 2 that the time inter- 
val between < t^5.6 refers to an inflationary period 
since the acceleration grows exponentially and the van 



der Waals fluid behaves like an inflaton with an equation 
of state p w = —p w . The time interval between 5.6<£<6 
is also another inflationary period since the acceleration 
is positive but it decreases and tends to zero whereas the 
energy density decays and the pressure is always nega- 
tive. For 6<i<6.4 the acceleration is negative and at- 
tains its minimum value whereas the pressure is positive 
and reaches its maximum value, hence this time interval 
can be interpreted as a matter dominated period. The 
next time interval between 6.4<i < oo refers also to a 
matter dominated period since the acceleration is always 
negative and the energy density and the pressure decay 
and tend to zero which corresponds to a dust dominated 
period. 

For 1/2 < w w < 2/3 the energy density decays with 
time, the pressure grows from a negative value to a max- 
imum positive value and decays at large values of time 
while there exists a period of acceleration followed by a 
period of deceleration. This interval for w w could sim- 
ulate an inflationary period (p w < 0), where the en- 
ergy density of the inflaton decays, followed by a matter 
(Pw 7^ 0) and dust (p w — ► 0) dominated periods. 

The energy density and the pressure have the same 
behavior as the former interval for 2/3 < w w < 3/4 but 
there exists only a period of deceleration and one can 
infer that this solution does not model properly the evo- 
lution of the Universe. 

For w w > 3/4 the pressure is always positive whereas 
the acceleration is always negative and one can conclude 
that it refers only to a matter (p w ^ 0) and dust (p w — > 0) 
dominated periods. 

Although the van der Waals fluid could describe the 
transition from an inflationary period to a matter domi- 
nated epoch it could not simulate the present accelerated 
period of the Universe. Moreover, in the present model 
we have not consider the irreversible processes concern- 
ing the energy transfer between the van der Waals fluid 
and the gravitational field. This will be the subject of 
the next section. 
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III. UNIVERSE AS A BINARY MIXTURE OF 
VAN DER WAALS FLUID AND DARK ENERGY 

If dissipative effects are taken into account and we 
consider the Universe as a binary mixture of a van der 
Waals fluid with dark energy, the energy-momentum ten- 
sor reads 

T» v = (p w +p d +p w + Pd + m)U»U»-( Pw + Pd + izj)g fw . 

(6) 

Above p d is the dark energy density, p d its pressure while 
w denotes the non-equilibrium pressure which is con- 
nected with the irreversible processes (see e.g. 0,0)- 

The evolution equation for the energy density of the 
mixture is now given by 

p w + pd + 3H(p w + p d + p w + Pd + tv) = 0, (7) 



while the Friedmann equation becomes 



U=Pu 



Pd- 



(8) 



The dark energy is supposed to interact only with itself 
and it is minimally coupled to the gravitational field, so 
that its energy density balance equation decouples from 
that of the van der Waals fluid, and we obtain from (JJJ 
two balance equations 



pd + 3H (p d + p d ) = 0, 



3H(p u 



Pu 



= -3Hm. 
(9) 

The term —3Hzu on the right-hand side of ©2 can be 
interpreted as the energy density production rate of the 
van der Waals fluid. It is connected with the energy 
density production rate of the gravitational field by (see 
e.g. i) 



p G + 3H(p G - p. u 



Pd) 



3Htu, 



(10) 



where pa is the energy density of the gravitational field. 
Hence, the non-equilibrium pressure is responsible for the 
irreversible transfer of the energy density of the gravita- 
tional field to the van der Waals fluid. 

Two models for the dark energy are considered here, 
namely the quintessence (p d — px) and the Chaplygin 
gas (p d — p c ). The equations of state for these two mod- 
els are given by 



Px = w x px, 
Pc = -A/pc, 



with wx < —1/3, 
with A = constant > 0. 



(11) 



For the motivation of the equations of state above one is 
referred to e.g. @>0] and the references therein. 

The equations of state (|llf> permit us to get from 
@i a relationship between the energy density of the 
quintessence and the cosmic scale factor as well as the 
corresponding one for the Chaplygin gas. In terms of 
dimensionless quantities these relationships read 



where Px/Pw an( l Pc/Pw refer to the initial amount (at 
t = by adjusting clocks) of the energy densities of the 
quintessence and of the Chaplygin gas with respect to 
the energy density of the van der Waals fluid, respec- 
tively and a is a dimensionless cosmic scale factor. The 
Chaplygin gas could interpolate a matter dominated Uni- 
verse (pressure-less fluid or dust) where p c oc 1/a 3 when 
ip/a 6 3> 1 and a cosmological constant dominated Uni- 
verse where p c — —p c when ?/>/a 6 <C 1. Hence, the pa- 
rameter ip - which is related to the integration constant 
of the differential equation ©1 - could be identified with 
the amount of the pressure-less fluid. 

The time evolution of the cosmic scale factor (0) for 
the mixture of a van der Waals fluid with quintessence 
becomes 
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Pa 



3{w x + l) 



= 0. 



(13) 



Equation (f 1 31) together with an equation for the non- 
equilibrium pressure lead to a complete determination of 
the time evolution of the cosmic scale factor for a mix- 
ture of a van der Waals fluid with quintessence. Within 
the framework of extended (causal or second-order) ther- 
modynamic theory the evolution equation for the non- 
equilibrium pressure in a linearized theory reads 



-3r)H 



or w 



(14) 



by assuming that the coefficient of bulk viscosity 77 and 
the characteristic time r are related to the energy density 
of the mixture p — p w + p d (here pd = px) by i] = ap 
with r = r\j p where a is a constant (see e.g. |(|). 

Once the cosmic scale factor is determined as a func- 
tion of time, the energy density of the quintessence can 
be calculated from (|12fh . while the energy density of the 
van der Waals fluid is obtained from 
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3(w x - 



(15) 



For the mixture of a van der Waals fluid with a Chap- 
lygin gas instead of l|13fl we have 



"4 
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= 0. 



(16) 



From the system of differential equations l|14J) and Ijl6(l it 
is possible to determine the time evolution of the cosmic 
scale factor for a mixture of a van der Waals fluid with 
a Chaplygin gas. Moreover, the energy density of the 
Chaplygin gas follows from ljT2|l 9 and the energy density 
of the van der Waals fluid can be determined from 



H 2 - 



1 



P° w VTTJj 



1 + 



(17) 



The system of differential equations i|13[l and l|14|) for 
the mixture of a van der Waals fluid and quintessence 
and the corresponding one (|16[) and (|14fl for the mixture 
of a van der Waals fluid and a Chaplygin gas have been 
solved by considering the initial conditions: a(0) = 1 for 
the cosmic scale factor, H(Q) = 1 for the Hubble param- 
eter and w(0) — for the non-equilibrium pressure. In 
order to have a complete determination of the time evolu- 
tion of the cosmic scale factor and of the energy densities 
there still remains much freedom, since each system of 
differential equations does depend on four parameters, 
namely: w w , p\j p^, wx, and a for the mixture with the 
quintessence and w w , Pc/p^, an d ot for the mixture 
with the Chaplygin gas. These parameters can be inter- 
preted as follows: (a) w w is connected with the inflaton, 
since the increase of w w leads to a less pronounced in- 
flationary period (see previous section); (b) p x /Pw and 
Pel Pw are relied to the initial amount of the energy den- 
sity of the quintessence and of the Chaplygin gas with 
respect to the energy density of the van der Waals fluid; 
(c) wx takes into account the strength of the negative 
pressure of the quintessence; (d) ib - as was previously 
commented - is related to the amount of the pressure-less 
fluid in the Chaplygin equation of state and (e) a is con- 
nected with the importance of the irreversible processes 
in the evolution of the Universe. 

In the figures 3 through 5 we have chosen w w = 0.6, 
p x /p° w = p° c /p° w = 0.03, w x = -0.9, i> = 3 and a = 0.13 
and we proceed to discuss the results. 

In the figure 3 it is plotted the "total" pressure as 
function of the time t for a Universe considered as a: 
(a) mixture of a van der Waals fluid with quintessence 
(Pw + Px + nj - straight line) ; (b) mixture of a van der 
Waals fluid with a Chaplygin gas (p w +p c + w- dashed 
line) and (c) one-component van der Waals fluid (p w +ru - 
dotted line) . In figure 4 it was adopted the same conven- 
tion for the lines in the plots of the acceleration field a as 
a function of time t. We infer from these figures that for 
the two models of the Universe in which the dark energy 
is present there exist: (a) an inflationary period with an 
exponential growth of the acceleration up to a maximum 
value where the " total" pressure remains with a constant 
negative value; (b) an interval where the acceleration de- 
creases from its maximum positive value to its maximum 
negative value corresponding to a growth of the "total" 




FIG. 3: "Total" pressure vs time for a mixture with quintes- 
sence (straight line), for a mixture with Chaplygin gas 
(dashed line) and for the one-component van der Waals fluid 
(dotted line). 




FIG. 4: Acceleration vs time for a mixture with quintessence 
(straight line) , for a mixture with Chaplygin gas (dashed line) 
and for the one-component van der Waals fluid (dotted line). 



pressure from its maximum negative value up to a max- 
imum positive value, and (c) a period where a growth of 
the acceleration takes place from its maximum negative 
value to a positive value and where the "total" pressure 
decays from its maximum positive value to a negative 
value. Almost the same conclusions can be drawn out 
here for the Universe modeled by the one-component van 
der Waals fluid, the difference being the lack of a negative 
" total" pressure and the corresponding accelerated epoch 
for the present period of the Universe, which is connected 
with the absence of the dark energy density. Hence, the 
mixtures of a van der Waals fluid with quintessence and 
with a Chaplygin gas can model the three periods of the 
Universe concerning its acceleration, beginning with an 
inflationary accelerated period, passing through a past 
decelerated epoch and leading back to a present acceler- 
ated phase, while the one-component van der Waals fluid 
can model only the two first periods, without an epoch 
of present acceleration. 

The energy densities of the van der Waals fluid p w , of 
quintessence px and Chaplygin gas p c together with the 
pressures of the quintessence px and of the Chaplygin 
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FIG. 5: Energy densities and pressures vs time for a mix- 
ture with quintessence (straight lines) and for a mixture with 
Chaplygin gas (dashed lines). 



gas p c are plotted in figure 5 as functions of time t. As in 
the figures 3 and 4 the straight lines refer to a mixture of 
the van der Waals fluid with the quintessence while the 
dashed lines correspond to a mixture with the Chaplygin 
gas. We conclude from figure 5 that the energy density 
of the van der Waals fluid decays more rapidly than the 
energy densities of the quintessence and of the Chaply- 
gin gas. This fact indicates that these two last energy 
densities outweigh the first one at later times during the 
Universe evolution, even if we consider a small amount of 
the energy density of the quintessence or of the Chaplygin 
gas with respect to the van der Waals fluid at the begin- 
ning (by adjusting clocks). We infer also from this figure 
that: (a) the energy density of the Chaplygin gas tends to 
a constant value at early times than the energy density of 
the quintessence does, and (b) the pressure of the Chap- 
lygin gas at later times is smaller than the pressure of the 
quintessence. These two observations are connected with 
the fact that the present acceleration of the mixture with 



the Chaplygin gas begins at early times than the corre- 
sponding one for the mixture with the quintessence (see 
also figure 4). 

Now we shall comment on the behavior of the solu- 
tions when the parameters for the mixture of the van der 
Waals fluid with the quintessence (w w , Px/Pwi w Xi ot) 
and for its mixture with the Chaplygin gas (w w , Pc/Pw, 
ip, a) are changed. The period of past deceleration be- 
gins at earlier times while the period of present accel- 
eration begins at later times when one of the following 
parameters changes: (a) the amount of the initial energy 
density of quintessence with respect to the van der Waals 
fluid Px/Pw decreases as well as the corresponding one 
for the Chaplygin gas p^/Pw: O 3 ) the strength of the neg- 
ative pressure of the quintessence wx reduces, and (c) 
the amount of the pressure-less fluid 4> m the Chaplygin 
equation of state increases. Hence, all these parameters 
are connected with a more pronounced predominance of 
the matter field in the decelerating phase of the Universe. 
The periods of past deceleration and the period of present 
acceleration begins at earlier times when w w increases or 
a decreases. This can be understood by recognizing that 
both parameters have influence on the early inflationary 
period, since the increase of w w leads to a less pronounced 
inflationary period while the decrease of a implies that 
at the beginning the non-equilibrium pressure has a less 
pronounced negative value. 

As a final comment let us write the equation of state 
of the Chaplygin gas as p c = w c p c . From figure 5 one can 
infer that for the times where the Chaplygin gas prevails 
over the van der Waals fluid we have that —1 < w c < 
—0.8. This last result together with the value wx = 
—0.9, adopted for the quintessence, are in agreement with 
the observational constraints for u^and wx presented by 
some authors (see e.g. the works 
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